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I. INTRODUCTION 

The semileptonic decays of B mesons are now well described with the aid of the 
heavy quark effective theory (HQET) [jl]-§l and the QCD perturbative corrections. 
On the other hand, the nonleptonic processes still suffer from lack of satisfactory 
theoretical description Q With this in mind, the former can be successfully employed 
in the determination of the parameters of the Standard Model (SM). Indeed they 
have been used to this end, recently yielding the Cabibbo-Kobayashi-Maskawa matrix 
element \Vcb\ [§]• Moments of lepton spectra can be used in the determination of a^, rrif, 
and rric P|-pT|.The particular processes involving the r lepton make this field yet more 
interesting as they are affected by the mass of the charged lepton, now comparable 
with the masses of the involved quarks mb and rric- 

The first order QCD correction to the differential decay width of the b quark 
has been found analytically |]T^ with the r polarization summed over. Now the 
polarization itself, which is no more suppressed by heavy quark masses, provides 
information on the SM parameters in its own right. It does not depend on the 
|Vca-m| elements, for instance, but can still yield the quark masses. This is all the 
more important due to the fact that the polarization turns out to be only weakly 
dependent on the coupling constant a^. It should be borne in mind that the first 
order QCD correction to the decay width itself is important, amounting to as much 
as 20% of the Born approximation. Thus the results of this paper render the r 
polarization especially applicable for use in evaluation of the quark masses. 

What we calculate here is the longitudinal polarization of the charged lepton. It 
may be added that the longitudinal polarization is transferred to the lepton system 
via the intermediating W~ boson whose longitudinal polarization in turn takes root 
in the Higgs mechanism, so that the process may shed light on physics beyond SM. 

We give the formulae for the decay of b quark into the charmed quark, r lepton and 
r-antineutrino in terms of the charged lepton energy and the squared four-momentum 
of the intermediating W~ boson, or, equivalently, the invariant mass of the r + z/ 
system. These formulae, combined with the ones for the unpolarized lepton case, 
readily give the polarization of the r lepton. This expression is then integrated to 
give r energy distribution. The correction to the polarization is shown together with 
the Born- approximated result. Then moments of energy distribution are evaluated 
for the polarized case. 

The sectioning of the paper goes as following. Sec. 2 is devoted to the kinematical 



'I'The accuracy of HQET for inclusive processes is still under debate, as the assumption 
of quark-hadron duality in the final state might introduce l/rric corrections not seen in 
the operator product expansion. Phenomenological analyses suggest that these dangerous 
terms are absent or small for two-quark processes like hadronic decays of r leptons and 
semileptonic decays of heavy quarks, see e.g. M and references therein. 



variables. In Sec. 3 the ideas behind the present calculation of the polarization are 
discussed. Then in Sec. 4 QCD corrections are briefly described. The flnal analytic 
result is given in Sec. 5 and in the following Sec. 6 the moments of energy distribution 
are given. 

II. KINEMATICS 

A. Kinematical variables 

In this section we deflne the kinematical variables used throughout the article 
as well as the constraints on those in both cases of a 3- and 4-body decay. The 
calculation is performed in the rest frame of the decaying b quark. In order to include 
the flrst-order QCD corrections to the decay, one must take into account both the 
3-body flnal state with a produced quark c, lepton r and an antineutrino Ur and the 
4-body state with an additional real gluon. The four-momenta of the particles are 
denoted as following: Q for the b quark, q for the c quark, r for the charged lepton, u 
for the corresponding antineutrino and G for the real gluon. All the particles assumed 
to be on-shell, their squared four-momenta equal their masses: 

,n2 2 2 2 2 2 2 /^2 n /t\ 

The four-vectors P=q+G and W=r + z/ characterize the quark-gluon system and the 
virtual intermediating W boson, respectively. The employed variables are scaled in 
the units of the decaying quark mass mb- 

_ml _ml _2E^ _W^ _ P^ 
m^ m^ mi, m^ mi 

Henceforth we scale all quantities so that ml = Q^ = 1. The charged lepton is 
described by the light-cone variables: 

-1 

Ti = -{x±^x^-4r)). (3) 

Thus, the system of the c quark and the real gluon is described by the following 
quantities: 

Po = lil-t + z) (4) 



Ps = ^Pi -z = ^[l + t' + z^-2{t + z + tz)]"^, (5) 

P^{z) = Po{z)±P,{z), (6) 

V ^ ^+^^) 1 ^+(^) (7^ 



where Po{z) and Psi^z) are the energy and the length of the momentum vector of the 
system in the b quark rest frame, yp{z) is the corresponding rapidity. Similarly for 
the virtual boson W: 

Woiz) = ^il + t-z), (8) 

Ws{z) = \lwi-t = ]^[l + e + z^-2{t + z + tz)Y'\ (9) 

W±{z)=W^{z)±W:,{z), (10) 

3;(.) = lln^^^ = ln^^^ (11) 

Kinematically, the three body decay is a special case of the four body one, with the 
four-momentum of the gluon set to zero, thus resulting in simply replacing z = p. 
The following variables are then useful: 

Po = Poip) = ^(1 - t + p) , P3 = Psip) = ^pI - P, (12) 

P± = P±{p) =Po±P3 , w± = W±{p) = I - p=p, (13) 

Y, = y,ip) = l\n^ , Y^ = yM = l^^—- (14) 

We also express the scalar products in terms of the variables used above, so in the 
units of the h quark mass one gets: 

Q-P=\{l + Z-t) ,T-V=]^{t-7^), 
1 1 

Q-v = -{I- z-x + t) , T- P = -{x-t-r]), 

Q.r = -X , iy-q= -{l-x-z + T]). (15) 

B. Kinematical boundaries 

The phase space is divided into two regions. The first of them, henceforth called 
A, is available for both three- and four-body decay, while the remaining part B cor- 
responds to pure four-body decay. Region A is defined as following: 

2y/ri <x <l + r] - p = Xm, (16) 

ti = r_(l - —^) < t < r+(l - —^) = t2. (17) 

1 — r_ 1 — T^ 

The additional region B of the phase space, where only 4-body decay is allowed, has 
the following boundaries: 

2y/rj<x<x^, r]<t<ti. (18) 



Conversely, if x should vary at a fixed value of t, the boundaries read: 

V<t<{l-y/pf, w. + ^<x<w+ + ^, (19) 

W- w^ 

for region A, and 

r?<t<v^(l-— ^), 2^<x<w^ + ^. (20) 

1 — ^Jrj w^ 

for region B. The upper limit of the mass squared of the c-quark-gluon system is in 
both regions given by 

Zmax = {l-r+){l-t/r+), (21) 

whereas the lower limit depends on the region: 

_ f P in Region A , . 

™" I (1 - r_)(l - t/r_) in Region B ^ ^ 



III. POLARIZATION AT TREE LEVEL 

In order to calculate the longitudinal polarization we must find the differential de- 
cay width for a given polarized final state of the r lepton. According to the definition 

r+-r" r^ 

P = TT, ;— = 1-2--, (23) 

r+ + r- r ' ^ ^ 

where F = r+ + F^. Once we know the width with the polarization summed over, 
we only need to find, for instance, the width for the negatively polarized r lepton. 
Before going to discuss the QCD corrections to the longitudinal polarization of the 
r lepton, we stop for a while to look at the tree level situation where the calculation 
is easy to follow. Once we choose the b quark rest frame and decide to look for 
the longitudinal polarization, we can express the lepton polarization four-vector s in 
terms of the four-momenta Q and r of the b quark and the r lepton, respectively: 

s = At + BQ (24) 

This is due to the fact that now only the temporal component of Q does not vanish, 
whereas the spatial parts of s and r are parallel. The coefficients A, B appearing 
in the formula above can be evaluated using the conditions defining the polarization 
four-vector s: 

s^ = -1 (25) 

s-T = 0. (26) 



Upon this one arrives at the following expressions: 

A^ = ±^^^, (27) 

B^ = ^^v^. (28) 

where the superscripts at A, B denote the polarization of the lepton. 

This observation combines with another one to make the whole calculation simpler. 
The total decay width dVo at the tree level for the unpolarized case reads: 

dVo = GlMl\VcKM?MZdn^{Q- q, r, iy)/7r' (29) 

where the matrix element amounts to 

MZir) = q-rQ-u. (30) 

With this kind of linear dependence on the four-momentum r, it is worth noting that 
the matrix element with the r polarization taken into account is, 

Mli = ^MZiK = T-ms) = ^{q ■K){Q. v\ (31) 

where m stands for the lepton's mass and we have introduced the four-vector K 

K = T-ms. (32) 



Applying now the representation (|2^ of the polarization s we readily obtain the 
following useful formula for the matrix element with the lepton polarized: 

-Mjs = T^^MZir) ± -^^—MZm (33) 

r 1 — r_ T< — T_ 



The first term on the right hand side of ( p3| ) can be calculated immediately once we 
know the result for the unpolarized case. Thus the problem reduces to performing this 
calculation again with the only difference amounting to replacing the four-momentum 
r of the lepton with that of the decaying quark, Q. The Born-approximated distri- 
bution can be written explicitly as 

r/r± 

-^ = 12ro/o^(x), (34) 

where 



Glm, 







F"%\Tr |2 



Vckm\. (35) 



1927r3 
The Born level function /o(a;) reads, for the unpolarized case. 



h{x) = ^C'r3 {C,[x^ - 3x(l + 77) + H + (3x - 6r7)(2 - x)] , (36) 



while the polarized cases are obtained using the function A/q: 

A/o(a;) = ^r|C'{C(3 - x - r/) + 3(x - 2)} (37) 

in the following way: 

/o^(^) = ^/o(a:)±A/o(x). (38) 

In the formulae above, 



rs = ^x^ - H C = 1 - , \ ■ (39) 

^ I — X + r] 

At the tree level, one can express the polarization integrated over the energy of the 
charged lepton as well: 

P = 1 - 1/18 [-2At]^ + {xl - 8r73/2)(3 + r/ - 3p) + 

12r/(a;^) - 3x^^/2 + 3(1 - r]f{p - p^/) + 

3(-l + 7])p{l - p/s^) [3(1 - r/)2 + p(3 + 57])] + 3^3/2 + 

3(x„ - 2^) (-1277 - V + 12r7p - 3p^ + 3rip^ + p^) - 

12r7p=^ In(sVp) - 18(27^2 - p^ - rj^p^) In [2^/{x^ + T3)] + 

18(l-r7Vln(2s2v^/[(l-77)(l-r7-T3)-p-r/p])} 

/ {T35/I2 + {2r,^ -p'- r]^p') ln[{xm + T-^/i^^)] 

-{l-r]^)pHn{2^p/[{l-r]){l-r] + T,)-p-r,p]}], (40) 

where 



s = l-^, T, = ^xl-4r], (41) 

^ = 1 - 7[(1 + v){v + P") + P(l + V^)] + r/^ + P^ + 12r/p. (42) 

Anticipating the subsequent discussion of the QCD corrections, let us already note 
that the specific linear dependence as featured in ( pO| , pTD goes back to the tensorial 
structure of the matrix element A^o.s: 

MZir) = C;:W\ (43) 

where £ and H stand for the leptonic and hadronic tensors, respectively. It is of 
course the leptonic tensor C where the linearity derives from: 

^t = T^^C;:{r) ± -^^Ci:{Q). (44) 

'^ r_|_ — r_ ^ r_|_ — r_ ^ 

It is not surprising then that the corrections to the hadronic tensor 7i will not affect 
this property. 



IV. CALCULATION OF QCD CORRECTIONS 

The Q CD-corrected differential rate for the 6 — > c + r_ + z/ reads: 

dT^ = dT^ + rfr^3 + (ir^4, (45) 

where 

dT^ = GlmllVcKMl'M^^^dn-siQ; q, r, z/)/7r5 (46) 

is the Born approximation, while 

rfr^3 = ^a,G'|m^|VfcKA/^-Mt3^7^3(g;g,r,^/)/7^6 (47) 

conies from the interference between the virtual gluon and Born amplitudes. Then, 

rfrj4 = '^asGlml\VcKM\^Mt^dn^iQ; q, G, r, u)/7i' (48) 

is due to the real gluon emission, G denoting the gluon four-momentum. Vckm is the 
Cabibbo-Kobayashi-Maskawa matrix element corresponding to the 6 to c or m quark 
weak transition. The Lorentz invariant n-body phase space is defined as 

dUniP^Pu ...,Pn) = S^'\P - EP^)n.^- (49) 

The superscript ± refers to the fact that now the polarization of the charged lep- 
ton is taken into account. In order to evaluate the appropriate rates it is conve- 
nient to take advantage of the decomposition (|4^ ) which led to the formula for the 
Born-approximated matrix element (0). As the QCD corrections influence only the 
hadronic tensor, the leptonic tensor and thus its linear dependence on r is left intact. 
This allows us to represent all the involved matrix elements in an analogous way 
and, in fact, the whole decay width is of the very same form: 

dr^ = T — rfr""(r) ± — ^ — rfr™(g). (so) 

r_|_ r_ r+ — r_ 

In Born approximation the contribution to the decay rate into the three-body final 
state is proportional to the expression 

_ 1 _ 1 

•^0,3 = 7^0 {x,t) = -q-KQ-u = 



'0,3 "40 v-^5 •-; — 2 
(1-p-x + t) 



Mx-t-ri)-ri{l+p-t)] (51) 



4(r+-r_; 
The three-body phase space is parametrized by Dalitz variables: 

dTZ'i{Q] q, T, u) = —dxdt. (52) 



The evaluation of the virtual gluon exchange matrix element yields: 

■^r3(^) = -[Hoq-rQ-u + H+pQ ■ uQ ■ r + H^q ■ uq ■ r 
+\p{H^ + H^)u ■ r + ip(/J+ -H_ + HL)\r-{Q-q- u)]{Q ■ v) 

-hi^\r.{Q-q-v)\{q.v) 

where 



(53) 



iJo = 4(1 - FpPoM) In Ag + (2poM)[^^2(l - 



p+w;^ 



-L%^{\ - — ) r,(y; + 1) + 2(in Vp + y;)(r^ + 1;)] 

+ [2p3i; + (1 - p - 2t) In ^]/t + 4, (54) 

i^± = ^[1 ± (1 - p)A]>;M ± ^ In v^, (55) 

i^L = ^(1 - In v/p) + ^ In ^ + |f,P3 + ^^. (56) 
and then the polarized case requires 

Ml, = -^^MZir) ^^MZiQ). (57) 

After renormalization, the virtual correction 7W^3 is ultraviolet convergent. How- 
ever, the infrared divergences are left. They are regularized by a small mass of gluon 
denoted as A^. In accordance with the Kinoshita-Lee-Nauenberg theorem,this diver- 
gence cancels out when the real emission is taken into account. The rate from real 
gluon emission is evaluated by integrating the expression 

where 

Bi{t) = q-r[Q-u{Q-G-l) + G-u-Q-uQ-G], (59) 
B^ir) = q- t[G ■ u - q- vQ ■ G + Q ■ u{q- G - Q ■ G -2q- Q)] 

+Q-u{Q-Tq-G-G-Tq-Q), (60) 

Bs{T) = Q-u{G-rq-G-pr-P). (61) 

Taking account of polarization amounts to substituting K for r in the coefficients 

^1,2,3- 

M^^ = ^+ MZir) — MY'M)- (62) 



The four-body phase space is decomposed as follows: 

dn^{Q; q, G, r, u) = dzdTZsiQ; P, r, z/)d7^2(P; q, G). (63) 

After employing the Dalitz parametrization of the three body phase space TZ^, and 
integration we arrive at an infrared-divergent expression. 

The method used in these calculations is the same as the one employed in the 



previous ones [0,|^,0. The infrared- divergent part is regularized by a small gluon 
mass Xg which enters into the expressions as In(AG')- When the three- and four-body 
contributions are added the divergent terms cancel out and then the limit Ac ^ 
is performed. This procedure yields well-defined double-differential distributions of 
lepton spectra as described below. 

V. ANALYTICAL RESULTS 

The following formula gives the differential rate of the decay b -^ tuX, X standing 
for a c quark or a pair of c and a gluon, once the lepton is taken to be negatively 
polarized: 



dV' { 12rn 



dxdt l2Tc 



Fo (a;, t) - ^FrA^^ t) for (x, t) in A 



2"-Fi7^(x,t) for (a;, t) in B 



(64) 



F{' differs according to which region (A or B) it belongs. Region A is available for 
the 3- and 4-body decay, while Region B with a gluon only. The following formulae 
are given for the negative polarization of the lepton, that is, we take 

A- = -^^^, (65) 

B~ = ^v^. (66) 

r I — r_ 



The factor Fq is defined in Eg. (|35|) , while 

F,-{x,t)= ^'^~P~''^^U n{x-t-v)-vil+p-t)] (67) 

T+-T. 

and 

F,;^(x, t) = Fo-$o + E D^^n + D^, (68) 

n=l 

F^^(x,t) = F,^o + E ^n^n + Di- (69) 



n=l 



The factor 12 in the formula (|6^) is introduced to meet the widely used pO|JT5|JTT 
convention for Fq{x) and Fq. The symbols present in (^) are defined as follows: 



$r 



2po 

P3 



[Lt2il 



t/r^ 



Ll2{l 



l-r+ 

) + Li2{l- 

P+ P+ 

l-r+. , l-t/r^ 



Ll2(l 



+Li2{w^) - Li2{w+) + 4Yp\n ^ 
Po. 



+4(l-^y;)ln(z 

P3 



Pj 4 m Zmax ; 

$1 = Li2{wJ) + Li2{w+) - Li2{T+) - Li2(t/r+) 

$3 = 2 In Vp, 

$4 = -ln(l-r+), 
$5 = ^ln(l-t/r+), 



P3 Z, 



max P \ , .1 / Zn 

+ 41n( — 



2po 

P3 



[Ll2{l 



- Ll2{l 

+Li2(l - 



l-r_ 



Lioil 



P+ 
1 -r_ 



P- 



mm P 
)+L22(l 

)-Ll2{l 

)+Li2{l 

-Li2{l- 
Lt2{r-) - 

^2 = - 



l-t/r+ 
P- 

P+ 

1 - t/r_ 
P+ 

-Lz2(t/r_) 
ln(l-r_) 



xI/3 = -ln(l-t/r_) 

^4 = -ln(l-r+) 

^5 = ^ln(l-t/r+) 



B. 



We introduce C1...C5 to simplify the formulae for D^ and D^: 
1 



(70) 
(71) 
(72) 

(73) 
(74) 
(75) 



(76) 
(77) 
(78) 

(79) 

(80) 

(81) 



Ci 



y/x'^ — At] 



p{—rix + Arjt — rjT^ — 4?7 + xr+ — tr^) — 2p rj + ri{xt + x 



10 



-At - 2t^ + 6 + 2r]) +T+{-r]x + r]t + 5r] - 2xt + x + t + f 



12) 



Co 



a 



p{v-t- 2r+) + {xt'^ - 2t^T+)/{2ri) + ri{-x - 2t + 2r+ + 30)/2 
+t(l +t-x) + T+{-2x + 2t + 6); 

[p{-lQrjx + 18?7t - 7r7r+ + xt + 7xt^ - lltr+ + 4r+) 



^3) 



v^aj^ — 4?] 



+p\-9r] + r+) + (2a;tV+ - xV/2 - t'rD/r] + r]\-T+ + 10) 
+xt{x -t- l)rj{Axt - XT+ - 2a; - a;^/2 + tr+ + At 



-Ar + llT^ 



12) 



r+(-6a;t - x + 23;^ - 3t + St^ 



C4 



p{2r]XT-^/t — Srjx + SrjT^/t + ISr^t — llr^r-i 



- 5)J ; {i 

- -A7^^T+/e 



■A) 



C. 






\/x^ — Ai] 

-xt + 5r+x - 7tr+) + p^{2riT+/t - 9r] - rfT+/t^) + xt(l + t - x) 
+ (-2a;tV+ + xH'^/2 + tV^)/^^ + rf{hT+/t^ - Qr^/t - 16/t + 6) 
+'q{2xT+/t + 2a;t - 3a;r+ - 4a; + a;V2 - 10r+/t + 3tr+ + 12t 
-4t^ + 3r+ + r^ + 4) + r+(-4xt - 3x + 2^^ + lit + 2t^)\ ; 
[p/{x - t - ri/t)]{-pri\+/t^ + riT+/t - r] - rfr^/t'' + rf/t + pr]T+/t 
-PV + prf/t)/'^ + [p/(l - a; + r7)](-r7t - ?7r+ + rf + tr+)/2 
+p(3r7r+/t + 9?7 - 9t - 3r+)/2 + {xt^/A - t^T+)/ri + ?7(3a; - 10r+/t 
-12t - 2r+ - 24 + 2r7)/4 + (3t + 5)r+/2 - a;t + 6t + 5tV2; 

1 



(85) 



^6) 

^7) 



{p?7(34 - 6a;r+/t^ - Ax/t^ + 3a;/t + 5a;t - 2a;r4 



A\/x^ — At] 

+3xyt^ + x^ + 8r+/t2 + 18r+/t - 2/t + 2tT^ 



2x 



38t - 14t^ + 20r4 



+pr7^(16 - 2x/t^ - 5T+/t^ - 2T+/t + 16/t - r+) + pr+(-4a;t - 4x + 6t 



+7t' 



11: 



P 



-^3 



r7(-10 + 6a;r+/t^ + 6a;/t^ + 4a;r+/t - 3a; - 3a; Vt^ 
-2a;Vt - 12r+/t2 - 18r+/t + 6/t + 18t - 6r+) + 7f{x/t^ + T+/t^ 
-T+/t - 2/t) + r+(2a; - 5t - 7) + +pri{-l{) - 2xT+/t^ - Ax/t^ - x/t 
+x^/t^ + 8r+/t2 + 6r+/t - 6/t) + prfT+/t^ + pr+ + p^'q{x/t^ - 2T+/t^ 
+2/t)] + r+(-4a;t + 2xt^ + 2x + 7t + t^ - 3t^ - 5) + ri{-lA + 2a;r+/t2 
+a;/t^ - 4a;r+/t - 2a;/t + 4a;t - 2a;t^ + 2a;r+ - a; - x^ /t^ + 2a;Vt - a;^ 

-2r+/t2 6r+/t - 10tr+ + 32t - 22t2 + At^ + 18r+) + 77^(28 + x/t^ 

-2x/t +x + 3r+/t2 - 5r+/t - 14/t + tT+ - lAt + r+)} ; (88) 

1 



2 /j.2 



^ : {p \r]{-6 - AxT+/t^ - ?,x/t^ - 2a;r+/t + x/t - 13a; + 2x^ /t 

\Jx^ — Ar] ^ L 

+x^/t + 6r+/t2 + 12r+/t - 2/t + 26t - 10r+) + 7f{-x/t^ - 2r+/t2 



11 






a 



A 



+2/t) +T+{6x + 6 - 12t)] + p^ri{-x/t^ + 2T+/t^ - 2/t) 

+p^ [r7(-10 + 2xT+/t^ + Sx/t^ - x'^/t^ - 6r+/t2 - 4r+/t + 4/t) 



Vr+/t 



+ rf{lA + x/t^ - x/t + 3r+/t^ - 2r+/t - 14/t - r_^ 



-1] T+/r - r+ + r7-(^i4 + x/r - x/t + 6T+/r - 2T+/t - , 
+r]{-U + 2xT+/t^ + x/t^ - 2a;r+/t - x/t + 2xt -2x- x^ jt^ + x^ jt 
-2r+/t^ - 8r+/t + 18t - 4^^ + 10r+) + r+(-2a;t + 2x + 2t + St^ - 5)} ; 

(89) 
-C^3 - C^\ (90) 

-C4 + C2; (91) 



4v^ 



= {[p/(l - X + r^)] \^{tr^ + t + r+) - 77^(1 + t + r+; 

+7^3 - tr+)] + [p/(x - t - r]lt)\{\ + p) \^{t - r+) + ^^^(-l + r+/t2+ 
r+/t - 1/t) +r73(l/t2 - r+/f'^)] + pr7(43 - a;r+/t - Axjt - 5a; + 2a;V^ 

+9r+/t + 13t + 5r+) + p^^{-^ - T+/t^ + 2r+/t + 1/t) + r+(9a;t + hx 
+4t - 6t^) + 12a;t + bxt^ - 2xH + p{-9xt - 3xr+ + 5tr+) + 
p^r7(-10 + 2a;/t - 3r+/t) + p^7f{-T+/t^ + 1/t) + (-6a;tV+ + xV 
+2tV+)/(2?7) + ?7(-35 - a;r+/t + 2a;/t - 2a;t + Axt+ + 26x 

-2a;Vt - 3a;V2 - 4r+/t 2tr+ - 34t - t^ - 22r+ - ^tD 

+if{-2 + 2x/t - X + 6r+/t + 2t)} ; 



(92) 



(93) 
(94) 
(95) 
(96) 
(97) 
(98) 
(99) 

One can perform the limit p — *> 0, which corresponds to the decay of the bottom 
quark to an up quark and leptons. The formulae, which are much simpler in this 
case, are presented in the same manner as the full results: 



1 


= Ci 


2 


= C2 - C3 


3 


= — C2 — C4 


f 


= C2 + C^ 


5^ 


= -c, + c. 




= c. 



with 



dT- 

dxdt 



F,-{x,t) 



12Tr 



Fq-Jx, t) - f^FiXx, t)\ for(x, t)inA 

for(x, t)inB 



12To^F,-^ix,t) 
(l-x-t) 



r+ 



T+{x-t-r]) -r]{l-t)] 



(100) 



(101) 



12 



and 



n=l 
_ _ 5 _ 

F,~b{x, t) = Fo^o + E ^n ^n + '^i- 



n=l 



where 



$0 = 2 



Li-) 



i^^2( ^, . )+^^2(t) 



2 



1-t ' ■ "^'^ 1/t-l 
+ ln2(l - r+)2 ln2(l - t) + ln2(l - t/r+) 
-21n(l -t)lnz^a^, 



TT 



•^1 = Y^ + Ll2{t) - L22(r+) - Lt2it/T+), 



$ 



203 



21n(l -t) 

$4 = $4, 
^5 = -^S, 



and 



^n = 2 



Li-} 



,T+-t,. 



Lin 



. l/r_-l , 

1-t ' ■ ""^' 1/t- 1 ■ 

l/r+ - 1 " 



Li^ 



,r+-t. 



-i:^2(- 



1/t-l 



ln(l - t) ln( 



1-t 



Z7-) 



--ln(^^)ln[(l-r+)(l-r_)] 
H\^^) In [(1 - t/r+)(l - t/r_)] , 



• 1 - t/T_ ■ 



^n = ^n, 



(n= 1...5). 



(102) 
(103) 



(104) 



(105) 

(106) 
(107) 



(108) 
(109) 



ct 



r+(x-2xt + t + t^ 



77(6 + a;t - a;r+ + a; + tr+ + 4t - 2t^ + 5r+) + 2?]^+ 



;iio) 



'-'2 



c.^ 



{xt^ - 2t^T+)/{2r]) + 77(30 -x-2t + 2r+)/2 - xt - 2a;r+ + 2tr+ 

+t + t2 + 6r+; (111) 

-^=L== [(4a:tV+ - xh' - 2t'rl)/{2r]) + +r]\lO - r+) 

r/(-24 + 8xt - 2xT+ - Ax - x^ + 2tr+ + 8t - 8t^ + 22r+ - 2r^)/2 



13 



(-4 



-6xtTjf 

1 



XT XT Xl~-L- \ X L ~r ^X 7~_L 



-3tr , + 5rr , 



5(-i^2: 



-4xtV+ + x^t^ + 2t^Tl)/{2ri) + ?7(4 + 2xT+/t + 2xt 



-3xT+ - 4x + xV2 - 10r+/t + 3tr+ + 12t 



At' 



3r , + rf 



+r/^(6 + 5r+/t^ - 6r+/t 16/t) - 4xtr+ + xt + xt^ • 

-x^t + 2xV+ + llir+ + 2tV+] ; 

(xt^ - 4tV+)/(4r7) + r/(-24 + 3x - 10r+/t - 12t - 2r+)/4 

+r/V2 - xt + 3tr+/2 + 6t + 5tV2 + 5r+/2; 



3xr4 






(113) 

(114) 
(115) 



P, 



203 






^ r7(-14 + 2xr+/t2 + x/t^ - 4xr+/t - 2x/t + 4xt - 2xt^ 

4vx"' — At] 

+2xr+ - X - xVi^ + 2xVi - x^ - 2T+/t' - Qr^/t - 10tr+ + 32t 

-22t^ + At^ + 18r+) + 7f{2S + x/t^ - 2x/t + x + 3r+/t2 - 5r+/t 

-14/t + tr+ ~14t + r+) + (-4xtr+ + 2xtV+ + 2xr+ + 7tr+ + tV+ 

-3tV+ -5r+)] 

— C2 — C3, 

C2 — C4, 



(116) 
(117) 
(118) 



:C5 + 



-6xrr+ + x^r + 2rr;)/(2r/) 



2"" ' 4Va;^ - 4?7 

+r7(-35 - xr+/t + 2x/t - 2xt + 4xr+ + 26x - 2x7^ - 3xV2 - 4r+/t 
-2tr+ - 34t - t^ - 22r+ - 3r^) + 77^ (-2 + 2x/t - x + 6r+/t + 2t) 
+ (9xtr+ + 12xt + Sxt^ + 5xr+ - 2x2^ + 4tr+ - 6tV+)] ; (119) 

On integration over t, one obtains tau energy distributions according to the formula 
I dV- „ , . 2a, 



/o"(^)-^/r(^)- (120) 

we obtain the corresponding formula for the 



12ro dx ■"' ' ' 3n 

where if we drop the superscript 
unpolarized case. 

The results are presented in Fig. |], where the polarization is plotted versus the 
r lepton energy. Both the Born and first order approximation are showed. In order 
to make the correction more explicit, we also present another diagram, where the 
function R{x) is drawn, defined as following: 

20- 

[l-Po(x)][l + -^i?(x)] (121) 



1 - P x) 



37r 



which gives the expression for R{x): 



R{x) 



fo{x) /o (x) 



(122) 
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FIG. 1. Polarization of r lepton in the Born approximation (dashed hne) and including 
the first order QCD correction (solid line) as functions of the scaled r energy x. The mass 
of the b quark taken at 4.75 GeV, c quark 1.35 GeV and the coupling constant a^ = 0.2 

thus its meaning is how the radiative corrections differ with respect to the state of 
polarization they act on. In ( |121| ) Pq denotes the zeroth order approximation to the 
polarization. The function fi{x)/fQ{x) has been presented in Fig. 0, while Fig. |^ 
shows the function R{x), so that one can immediately see how small the correction 
is in comparison to the correction to the energy distribution. 

Integrating over the charged lepton energy one can obtain its total polarization 
as well as the corrections to which it is subject. If we take rrif, = 4.75 GeV as the 
central value for the decaying quark mass and rUb = 4:A and nif, = 5.2 for the limits, 
we arrive at the following (the mass difference rrif, — rric = 3.4 GeV everywhere): 




(1 - Po) 1 + -^Rs + -^RL + —A (123) 



with 



-To — -U. M88+0.0109 

Rs = — 0.016^q;oi7 

Dl _ Q 421+0.027 

-'^np ~ '-'•^^-'--0.025 
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FIG. 2. The ratio f^ (x)//q (x) representing the radiative correction for the negatively 
polarized state as dependent on the scaled r lepton energy x. 

7~)2 o r)Q— 0.22 

^np — ^-^"+0.16 



VI. MOMENTS OF r ENERGY DISTRIBUTION 

The moments of r energy distribution, which are useful sources of information on 
the physical parameters regarding the discussed decay can be evaluated according to 
the formula: 

Mit = / E^-T^dE^, (124) 

'^ Efrtin (111/ J- 

'■" = w <^2^) 

where Emm and E^ax are the lower and upper limits for r energy and M„ include 
both perturbative and nonperturbative QCD corrections to r energy spectrum. The 
superscripts denote the polarization states. Since one obviously has 

M^ = M^ + M- (126) 

where M„ stands for the unpolarized momenta, we only give the values of the momenta 



for the negative polarization case. The unpolarized distributions were given in [12 
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FIG. 3. The QCD-correction function R(x) for the pole mass values of the b quark 
taken to be 4.4 GeV (dashed), 4.75 GeV (solid) and 5.2 GeV (dash-dotted) as dependent 
on the scaled r lepton energy x. 



The nonperturbative corrections to the charged lepton spectrum from semileptonic B 
decay have been derived in the HQET framework up to order of 1/ml |[T5|-[T7|. The 
corrected heavy lepton energy spectrum can be written in the following way: 



1 dT 
UTodx 



Sir m 



2 J np V 






m 



2 J np 



X) 



(127) 



where Ai and A2 are the HQET parameters corresponding to the b quark kinetic 
energy and the energy of interaction of the b quark magnetic moment with the chro- 
momagnetic field produced by the light quark in the meson B. The functions f^h'^^ 
can be easily extracted from the formula (2.11) in [^].The formula (|127|) looks identi- 
cally if one considers definite polarization state of the final r lepton. The appropriate 
calculation within the HQET scheme has also been performed |jl5[, see formula (2.12) 
therein. 



Following ref. |Tl]] we expand the ratios 



,(0)- 



on '^ 2 n 2 n 

in vat vat 



128) 
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where r^°) is the lowest approximation of ■ 



,(0)- 






,/. 



l+ri p 

2v^ 



/o {x)x^dx 



J 2^1 Jo {x)dx 



(129) 



Each of the 5^*^ is expressed by integrals of the corresponding correction function 
/'•*•' (x) and the tree level term fo{x) 



5«- 



/2^J|-— ^/W-(x)a:"rfx /2^ ;|-— ^/W-(x)rfx 

i2,A7 /o (a^)"a^ 



(130) 



J2^ /O [Xjx'^dx j2^ 

where the index i denotes any of the three kinds of corrections discussed above. The 
coefficients 6!^^ depend only on the two ratios of the charged lepton and the c quark 



to the mass rrib. Following ref. [|12| we employ the functional dependence of the form 

5«-(m,,m„m.)=5«-f^,^V (131) 



\ rrir rUb J 

The quark masses are not known precisely so we have calculated the coefficients 
in a reasonable range of the parameters, that is, 4.4 GeV< rUb < 5.2 GeV and 
0.25 < iTLc/mh < 0.35 and then fitted to them functions of the following form: 

5(p, q) = a + b{p- po) + c{q - go) + d{p - pof 

= +eiP-Po)iq-qo) + fiq~qo)^ (132) 

where p = mf,/mr,po = 4.75 GeV/1.777 GeV= 2.6730, g = mc/mbjqo = 0.28 and 
the polynomial coefficients can be fitted for each of the (J^*-* separately with a relative 
error of less than 2%. Our choice of the central values reflects the realistic masses of 
quarks: rub = 4.75 GeV and rric = 1.35 GeV,for which 5^*^ = a^\ 

To bring out the difference in the extent to which the corrections affect the two 
different polarization states it is useful to compare these coefficients with the ones 
obtained with the polarization summed over. This can be done along the lines sug- 
gested by the treatment of the polarization itself, see Eq.( |121| ). The corresponding 
expansion takes the form 



r 

n 



-(0) 



rl°) 



l-|^(.^)--.^)) + ^( 



(1)- 



A 



^^') + -M'- - ^^') 



mt 



(133) 



With these, one can readily find the actual relative correction of each kind, assuming 
reasonable values of as,Xi andA2. Here we take as = 0.2, 0.15 GeV^ < — Ai < 
0.60 GeV^, A2 = 0.12 GeV^, keep the b quark mass fixed at 4.75 GeV and the mass 
of the c quark equal to 1.35 GeV. The corrections then read for n = l{n = 5) 





Correction to 


Type 


r- /A0)- 


rn/r^r^'^ 


perturbative 


-0.00084(-0.0048) 


-0.0009(-0.0052) 


kinetic energy 


0.008 ±0.005(0.06 ±0.04) 


0.008 ±0.005(0.06 ±0.04) 


chromomagnetic 


-0.0097(-0.053) 


-0.0092(-0.0511) 
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VII. SUMMARY 

The first order perturbative QCD correction to the polarization of the charged 
lepton in semileptonic B decays has been found analytically. It is expressed by the 
charged lepton energy and the invariant mass of the lepton system. The polarization 
correction has turned out to be very small as it does not exceed 1% taking common 
values for the parameters occuring in the formulae. This makes the polarization a 
very useful quantity in determining the quark masses. The moments of r lepton 
energy distribution have been evaluated for the case of a polarized lepton and the 
correction has again been found to be little different from the one in the case of an 
unpolarized lepton. The nonperturbative HQET corrections to the moments have 
also been calculated using the formulae from ||15| |. 
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